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ABSTRACT 

We compare the luminosity functions for red galaxies lying on the restframe (U — V) color-magnitude 
sequence in a homogeneous sample of ten X-ray luminous clusters from the MACS survey at z ~ 0.5 
to a similarly selected X-ray cluster sample at z ~ 0.1. We exploit deep Hubble Space Telescope 
ACS imaging in the F555W and F814W passbands of the central 1.2-Mpc diameter regions of the 
distant clusters to measure precise colors for the galaxies in these regions and statistically correct for 
contamination by field galaxies using observations of blank fields. We apply an identical analysis to 
ground-based photometry of the z ~ 0.1 sample. This comparison demonstrates that the number of 
faint, My ~ —19, red galaxies relative to the bright population seen in the central regions of massive 
clusters has roughly doubled over the 4 Gyrs between z ~ 0.5 and z ~ 0.1. We quantify this difference 
by measuring the dwarf to giant ratio on the red sequence which increases by a factor of at least 
2.2 ± 0.4 since z ~ 0.5. This is consistent with the idea that many faint, blue star-forming galaxies in 
high density environments are transforming onto the red sequence in the last half of the Hubble time. 
Subject headings: Galaxies: clusters - galaxies: luminosity function - galaxies: evolution 
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1. INTRODUCTION 

Recent studies have precisely quantified the variation 
in the photometrically classified galaxy po pulation as 
a fun ction of environment at low redshifts (lHogg et al.l 
2004). These studies separate galaxies into red, pas- 
sive or blue, star-forming systems and find that the 
proportion of the latter decr eases in higher den sity re- 
gions in the local Universe (|Baldrv et al.l [20061. The 
identification of the physical process responsible for this 
trend is still contentious, in part because it is likely 
that a number of processes contribute in different en- 
vironments, at different epochs and acting on galax- 
ies of different luminosities/masses. The presence of a 
range of potential evolutionary pathways linking star 
forming and passive galaxies may be reflected in the 
diversity of star formation histories derived for pas- 
sive early type galaxies. While the formation of lu- 
minous early type galaxies in clusters has been inter- 
preted in terms of a narrow range in star formation his- 



Aragon-Salamanca et al.lll993t 
Stanford et al.lll998ft . there is 



tories (Bower et al.111992) : 
Ivan Dokkum et alj 11998c 
evidence of much more variety in lo wer luminosity 
systems (see iFerguson fc Binggelil 11994 for a review) . 
Several lines o f evide nce illustrate this, for example 
iPoggianti etall (|2001h find a broad range in ages but 
a slight decre ase in age for fainte r dwarf galaxies in the 
Coma cluster. iSmail et al.l (|2001f ) reached a similar con- 
clusion for the luminosity weighted ages of low luminosity 
early type galaxies in the z = 0.18 cluster A 2218. An 
even wider variety is found in much lower luminosity sys- 
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terns, iGrebell (| 19991 ) finds that star formation timescales 
and ages vary significantly for the local group dwarf ellip- 
ticals. Indeed, this same effect may underlie the varying 
morphological mix seen in the passive galaxy populations 
in clusters, where there is a claim of a deficit of early 
type disk galaxies (pred ominantly SO) in distant clus- 
ters ([Dressier et al.lll997l ). Thus it appears that passive 
galaxies may be formed via a number of different pro- 
cesses and that the mixed nature of the population may 
be most easily discerned at the lowest luminosities. 

These different pathways may also result in a chang- 
ing passive galaxy population in clusters at different red- 
shifts. The most fundamental measure of the transfor- 
mational processes forming passive galaxies in high den- 
sity regions is to look at t he build up of th e luminosity 
function of this population. iDe Lucia et al.1 (|2004l ) there- 
fore investigated this scenario by measuring the color- 
magnitude relation in four optically selected z ~ 0.75 
clusters and comparing these to the nearby Coma clus- 
ter. They find that the high redshift clusters exhibit 
a deficiency in lo w luminosity red gala xies compared to 
Coma. Similarly, iKodama et al.l (|2004h found a deficit of 
red sequence galaxies when looking at z ~ 1 high den- 
sity regions in the Subaru/ XMM-Newton Deep Survey. 
This leads them to conclude that many faint red galax- 
ies in clusters have onl y moved onto the red sequence 
since z ~ 0.75 (see al solDe Lucia et all 12007). This re- 
sult is controversial as lAndreonl (|200rf " claims there is no 
evidence for this red sequence build up when comparing 
the luminosity function of a cluster at z = 0.831 to those 
of local clusters. 

In this paper we aim to test these results by compar- 
ing the evolution in the luminosity function of galaxies 
on the red sequence in two well-defined samples of X- 
ray selected clusters at z ~ 0.5 and z ~ 0.1. We em- 
ploy Hubble Space Telescope (HST) optical imaging of 
ten z ~ 0.5 X-ray luminous clusters and compare these 
to a similar sample at z ~ 0.1 to examine evolution in 
the faint end of the red sequence luminosity function. 
We adopt a ACDM cosmology (fi M = 0.3, Sl Vac = 0.7, 
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Hq = 70kms~ 1 Mpc -1 ) in which the lookback times to 
z = 0.13 and 0.54 are 1.6 and 5.3 Gyrs and the angular- 
scales are such that 1" corresponds to 2.3 and 6.4 kpc re- 
spectively. An AB magnitude system is used throughout 
this paper. 

2. OBSERVATIONS AND REDUCTION 

Our analysis employs restframe optical imaging of two 
cluster samples at z ~ 0.1 and z ~ 0.5. The z ~ 0.5 
sample are selected from th e MAssive Cluster Survey 
(MACS, lEbeling et afll200l . MACS is a survey of dis- 
tant X-ray luminous, (Lx > 10 44 ergs _1 ) and therefore 
massive, galaxy clusters selected from the ROSAT All- 
Sky Survey. The ten clusters in this sample, along with 
a further two for which archival HST observations ex- 
ist, are the twelve most distant clusters, z ~ 0.5-0.7, 
from MACS and all have X-ray luminosities of L(0.1- 
2.4 KeV) > 15 x lO^ergs" 1 (Table 1). These ten clus- 
ters were imaged with the Advanced Camera for Sur- 
veys (ACS) Wide Field Channel (WFC) on HST dur- 
ing Cycle 12 (GO#9722). Each cluster was imaged for 
two orbits (4.5 ks) through both the F555W (V 555 ) and 
F814W (I814) filters. These data were retrieved from the 
STScI archive and reduced using standard STScI soft- 
ware (multidrizzle v2.7) to provide final images with 
0.05" sampling and good cosmetic properties. 

We extracted the photometry from the AC S images us- 
ing SExtractor (|Bertin fc Arnoutsl[l996f ) run in dual- 
image mode so that photometric information from the 
Vg55-band was extracted for all sources detected on the 
J8i4-band image and ensuring that we have precise aper- 
ture photometry for even crowded sources. Throughout 
this paper we use 1.4-" apertures (9 kpc diameter) to cal- 
culate the (V555 — Igu) color and the 'Best' magnitude is 
used for the total isi4-band magnitude (these correspond 
closely to restframe (U — V) colors and V absolute mag- 
nitudes). 

The low-redshift, z ~ 0.1, comparison sample for 
our analysis comes from t he Las Campanas/AA T Rich 
Cluster Survey (LARCS. IPTmbblet et all l200l 2006). 
This survey obtained panoramic, ground-based B- and 
i?-band imaging of ten X-ray luminous (Lx > 5 x 
10 44 ergs _1 ) clusters at z = 0.08-0.15 selected from the 
ROSAT All-Sky Survey (Table SJ. The observations and 
their reduction and a nalysis are described in detail in 
IPimbblet et~aT] (|200lD . Here we use 4-" (B - R) colors 
(corresponding to restframe (U — V) colors in 9 kpc aper- 
tures) and total magnitudes derived from these restframe 
~ V^-band selected galaxy catalogues in our analysis. 
These data are thus well-matched to the observations of 
the distant sample: the absolute T^-band surface bright- 
ness limits are \iy — —15.4 and —15.8 mag. arcsec -2 for 
LARCS and MACS respectively, with spatial resolution 
of 2.8 and 1.0 kpc. In our analysis we only consider the 
inner parts of each cluster, within a radius of 600 kpc of 
the cluster center as identified from the X-ray emission 
(usually corresponding to the position of the brightest 
cluster galaxy) as this provides uniform coverage across 
both the LARCS and MACS datasets. 

Both the MACS and LARCS cluster samples were se- 
lected from the same X-ray survey and the luminosi- 
ties for the clusters are sufficiently bright that they 
should correspond to some of the most extreme envi- 
ronments at their respective epochs. The median X- 



ray luminosities of the high- and low-redshift samples 
are 17.0 and 7.3 x 10 44 ergss _1 respectively, correspond- 
ing to a difference of less than a factor of two in the 
typical total mass dReiprich fc Boehringer|[2002T) . How- 
ever, an important issue to address is that the mass of 
the z ~ 0.5 progenitors of the LARCS clusters would 
be eve n lower than th e MACS clusters. Using the re- 
sults of lTormenl (|1998f) we see that the progenitors of the 
LARCS clusters at z ~ 0.5 would be ~ 3.5 x less massive 
than the MACS sample (with c orresponding X-ray lu- 
minos ities of ~ 2 x 10 44 ergss~ 1 . lReiprich fc Boehringerl 
(2002)). There is no evidence for strong variations in 
the galaxy luminosity function between clusters span- 
ning such a relatively m odest differe nce in typical mass 
(|de Propris et al.lll999fh In addition. IWake et alj (|2005f ) 
see no variation in the blue galaxy fraction in clusters 
covering a similar range in mass. Combining these two 
results, we therefore expect that any differences between 
the galaxy populations in these two samples will primar- 
ily reflect evolutionary differences between z ~ 0.5 and 
z~0.1. 

Finally, to better define the evolutionary trends we 
are searching for, we also include similar observations 
of seven additional clusters in our analysis of the dwarf- 
giant ratio in ^3.31 We include a low redshift point from 
U - and V^-band observations of the Coma cluster from a 
dataset of known members (jGodwin et al.l [19831 ). This 
is well matched to our main MACS and LARCS sam- 
ples. Next, we include four additional clusters which are 
at redshifts intermediate between the LARCS and MACS 
samples and a further two clusters at redshift beyond the 
MACS sample (see Table[T]). These clusters all have deep 
archival observations with HST in the F555W/F606W 
and F814W filters (from GO projects 9033, 9290, 9722, 
10325, 10491, 10872 and 10875) and we obtain the data 
from the HST archive. We note that the HST ACS ob- 
servations of the four intermediate clusters do not reach 
the 600-kpc radius adopted for our analysis and so we 
have corrected the dwarf-giant ratio using the average ob- 
served radial trend in this ratio for the composite MACS 
sample. The dwarf-giant ratio is observed to increase 
in value from the cluster center outwards and flatten at 
larger radii. For the most extreme case, cluster Abell 
1703 at z = 0.258, this correction gives an enhancement 
of 15% in dwarf-giant ratio. In addition we caution that 
the K-corrections for the additional HST clusters are not 
as well- matched to restframe (U — V) as those for the 
LARCS or MACS samples and so there may be system- 
atic uncertainties related to these data points. 

3. ANALYSIS AND DISCUSSION 

We show in Fig. Q] the color-magnitude diagrams 
for the individual MACS z ~ 0.5 clusters (similar 
plots for the individ ual LARCS clusters are given in 
Pimb blet et al.ll2002f ). These display strong linear fea- 
tures in the distributions of the brighter and redder 
galaxies in the fields. These are the color-magnitude 
relations for the passive, ear l y-type cluster me mbers 
(jSandage fc Visvanathanl 119781 : Bower et alj|1992f ). For 
the brighter galaxies in the clusters this sequence is 
thought to represent variations in the metallicity and to a 
lesser extent age of the stel l ar populations in the galaxies 
(|Kodama fc Arimotdll997l : iTerlevich et al.ll2001l ). 

As our MACS and LARCS cluster samples are ho- 



An Increase in the faint red galaxy population in massive clusters since z~0.5 



3 



mogeneously selected and observed and cover a small 
range in rcdshift (az/z ~ 0.15), we are able to com- 
bine them to reduce the influence of variations in field 
contamination and to improve the statistics of our anal- 
ysis. We show in the two lower panels of Fig. [T] the 
combined color-magnitude relations for the MACS and 
LARCS samples. The combined color-magnitude plots 
were created by correcting the data to the median red- 
shifts of the MACS and LARCS samples, z = 0.54 and 
z = 0.13 respectively, u sing the K-corrections from a 
iBruzual fc Chariot! (|2003f ) solar metallicity, simple stellar 
population, Zf = 5 model and the appropriate distance 
modulae. We also define limits to the color-magnitude 
relation in these combined samples to allow us to quanti- 
tatively compare them. We define the limits on the basis 
of an error weighted two-parameter fit to the slope of the 
combined MACS red sequence with a stripe width of 0.75 
mags to include the observed scatter. The corresponding 
boundaries for the LARCS red sequence are then deter- 
mined by K-correcting and c olor converting the limits 
from the MACS sample using IBruzual fe Chariot! (12 003) 
and the formulae given in Skiff (120031) and lNatali et al.l 
(1994). A correction for the observed change in red se- 
quence slope between the two epochs is also included by 
using the gradient found from error weighted fitting to 
the combined LARCS red sequence. We plot on Fig. [1] 
the corresponding color boundaries for the two samples. 
It is these red sequences that are used below to estimate 
the combined red sequence luminosity functions and the 
relevant field correction. 

3.1. Field Correction 

The removal of field galaxies from our samples is crucial 
to provide a clean measurement of the cluster luminos- 
ity function. For the MACS HST sample we used blank 
fields from The G reat Observatories O rigins Deep Survey 
North (GOODS. iGiavalisco et afll200l . These provide 
estimates of the correction for field contamination in each 
bin in our color-selected luminosity functions. The vari- 
ation in this correction between independent 16 arcmin 2 
subregions of the 112 arcmin 2 blank field is included in 
the luminosity function errors. The correction for field 
contamination for the ground-based observations of the 
LARCS sample is calculated in a similar way. Here, 
however, we make use of the fact that the panoramic, 
2-degree diameter, imaging of these clusters extends into 
the field population surrounding the clusters and so we 
can use the outskirts of the images to determine the field 
contamina tion. These estimates have been shown to be 
robust bv iPimbblet et alj (|2002t ). A gain we determine 
the reliability of our field corrections by dividing the to- 
tal 4600 arcmin 2 field region into 290 arcmin 2 subregions 
and determining the scatter in these independent areas. 
These field subregions are on scales comparable to the 
regions of the MACS and LARCS clusters we analyse. 
The typical 1-sigma errors for the field correction are in 
the region of 10 percent for MACS and 20 percent for 
LARCS. We propagate these uncertainties through our 
entire analysis. 

An additional source of contamination comes from 
higher redshift galaxies that have been gravitationally 
lensed by our clusters. This would effectively alter the 
field contamination in our sample, either increasing or 
decreasing it depending upon the slope of the number 



counts (|Tavlor et al.l [1998) . To estimate the scale of 
this effect we compare the number counts of galaxies 
on the color-magnitude plane which are just redward of 
the cluster's color-magnitude sequence (and hence should 
be at higher redshifts) to similar regions of the color- 
magnitude plane for the blank fields. This provides an 
estimate of the potential enhancement in the surface den- 
sity of red galaxies in the cluster centres due to lensing 
of ~0.1 percent. Although only a very small effect, we 
include this factor in the field correction. 

3.2. Luminosity Function 

The field-corrected luminosity functions for the two 
composite cluster samples are shown in Fig. [2j Lu- 
minosity functions are traditionally fitted with a single 
Schechter function (|Schechterlll976h . Recent papers on 
clusters, however, have instead fitted a Gaussian to the 
bright end of the luminosity function and a Schechter 
function to t he faint end as these give a better fit to 
obser vations (|Thompson fc Gregorvlll99l : iDahlen et all 
l200l . We compare the single Schechter to the Gaus- 
sian+Schechter parameterisation of the luminosity func- 
tion for galaxies on the red sequences in the LARCS and 
MACS samples. To avoid incompleteness we only con- 
sider the luminosity function down to the K-corrected 
5-cr limit of the highest redshift cluster in the shallower 
LARCS data (R = 20.92 from A 3888 corresponding to 
M v = -17.75). Table H contains the best-fit Schechter 
and Gaussian parameters converted to absolute y-band 
magnitudes using the method described above and the 
reduced x 2 f° r these fits. The errors quoted here are 
estimated using a bootstrap method. 

When fitting combined Gaussian+Schechter functions 
we fix the mean magnitude and width of the Gaus- 
sian components (Table H]), within their errors, so as 
to constrain the evolution of the bright-end Gaussian 
between z ~ 0.5 and z ~ 0.1. We do this as the lu- 
minosity evolution of galaxies in the bright end of the 
LF is well-constrained fr om fundamental plane studies 
(|van Dokkum et aTJll998[ ) and so we can focus on changes 
in the faint-end. The passive evolution of the luminosity 
between the two epochs is fi xed as 0.3 magnitudes from 
Ivan Dokkum fc Franxl (|200lD . 

We find that both a Gaussian+Schechter or a Schechter 
function give acceptable fits to our two samples. Both 
parametric forms also demonstrate the same difference 
between the two samples: an increase in the number 
of faint red galaxies compared to the brighter red pop- 
ulation at lower redshifts. For the Schechter fits this 
is shown by the steeper faint end slope (a) in the 
LARCS clusters than MACS clusters, -1.11±0.06 versus 
—0.91 ± 0.02 respectively, a difference of approximately 
3.2a. We find no evidence for variations of the form of the 
luminosity function looking at either the bluest or red- 
dest halves of the color-magnitude sequence or between 
different clusters when ranked by richness. 

For the Gaussian+Schechter fits, the change in the lu- 
minosity function is shown in part by the relative nor- 
malisation of the faint Schechter and bright Gaussian 
components, </>*/Amp, which decreases from 1.94 + 0.58 
in LARCS to 1.26 + 0.67 for the MACS sample. However, 
the covariance between the parameters in the two compo- 
nents functional fits makes such a comparison complex to 
interpret and so we turn to another, shape-independent, 
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parameterisation of the relative numbers of faint and lu- 
minous galaxies: the dwarf-giant ratio. 

3.3. Dwarf-Giant Ratio 



The results from £)3.2I are difficult to interpret in part 
because the form of the luminosity function is complex 
and its evolution is uncertain. A far simpler approach to 
quantify the relative evolution of the numbers of bright 
and faint galaxies is to use the ratio of the number of 
dwarfs to giants along the red sequence: the Dwarf-Giant 
Ratio (DGR) which provides a single number to param- 
eterize the distribution of galaxy luminosities within a 
population. The variation of this quantity (or its inverse 
GDR) with distance from cluster center, density and 
cluster richness have been studied by a number of worker s 
(jDriver et al.lll998t iDahlen et aill2004l ; iGoto et al.ll200l . 
Therefore to provide a shape-independent comparison of 
the red galaxy populations in the MACS and LARCS 
clusters we have measured the DGR. 

The boundary between giants and dwarfs is arbi- 
trary and is usually defined as the magnitude where 
the faint-end Schechter fu nction begins to d ominate over 
the bright-end Gaussian ( Goto et al.l l2005h . Looking at 
the distributions in Figt2l we therefore choose an ab- 
solute magnitude of My ~ —19.9 at z = 0.13 as our 
dividing point. This absolute magnitude brightens to 
My ~ —20.2 at z = 0.54 as we take into account the 
passive evolution models of Ivan Dokkum fc Franxl (|2001f ) 
(we have confirmed that our results are not sensitive to 
applying this factor). To ensure that our measurements 
of the DGR in the different samples are not effected by 
incompleteness we adopt the same faint-end limit as for 
fitting the luminosity functions (My ~ —17.75) and as 
in §3. 21 we only consider galaxies within a 600 kpc radius 
of the center of each cluster. The limits chosen for our 
PGR analysis are comparable to those of lDe Lucia et all 
(|200l . 

Fig. [2] shows the variation of the DGR on the red se- 
quence (RDGR) with redshift. The RDGR increases 
with cosmic time and we attribute this to an increas- 
ing number of dwarfs on the red sequence. We deter- 
mine weighted mean RDGRs for the MACS sample of 
1.33 ± 0.06 and 2.93 ± 0.45 for LARCS, a difference of 
3.7<7. This corresponds to an increase in RDGR of a 
factor of 2.2 ± 0.4 from z = 0.54 toz = 0.13, or a look- 
back interval of 4 Gyrs. An alternative way to look at 
this evolution is the variation in the proportion of inte- 
grated red light at the faint-end of the color-magnitude 
relation, which increases by a factor of 1.46 ± 0.14 from 
z = 0.54 to z = 0.13. This means that the stellar mass 
in the passive dwarf population (My < —20) now al- 
most equals that in luminous cluster galaxies, whereas 
at z ~ 0.5 the giants dominated the total F-band lu- 
minosity from galaxies on the color-magnitude relation. 
The errors shown are a combination of the Poisson un- 
certainty and the field correction error. 

We note that the LARCS RDGR errors are larger than 
those for MACS mainly due to the fact that at z ~ 0.1 
the field galaxies and the faint end of the red sequence 
occupy the same region of (B — R)-R color-magnitude 
space increasing the field correction error at the faint end 
(see Fig. [IJ. In contrast the faint red sequence in the 
MACS sample at z ~ 0.5 is much easier to distinguish 
from the field using V and / bands. 



To parameterise the evolution in the RDGR we fit a 
(1 + z)' 13 power-law to the LARCS and MACS samples 
in Fig. O We see that this provides a good descrip- 
tion of the evolution for f3 = 2.5 ± 0.5, with all of the 
clusters consistent with the fit. This confirms that the 
luminosity function of the red sequence in the central re- 
gions of massive clusters appears uniform with no clear 
evidence from our study of real cluster-to-cluster varia- 
tions although the errors on individual clusters are large. 
We also find that the six additional clusters and Coma 
(RDGR= 2.8 ± 1) follow th e same trend defined b y the 
MACS and LARCS samples. iDe Lucia etai] (|2004| ). who 
use a similar definition of DGR, found a value of 1.23 for 
clusters at z = 0.75 and a value of 2.9 for Coma which 
are in good agreement with the trend we observe. The 
trend is also in quali tative agreement with the work of 
IKodama et ail (|2004h . 

4. CONCLUSIONS 

Our analysis of the red galaxy populations in X-ray 
luminous clusters shows clear differences in the form of 
the luminosity function over the redshift range z = 0.1 - 
0.5. These changes reflect an increase in the proportion 
of dwarf to giant galaxies in the population since z ~ 0.5 
which we attribute to an increase in the number of dwarfs 
on the red sequence. We quantify this evolution using the 
shape independent estimate of the relative evolution of 
the faint end of the luminosity function, the red sequence 
Dwarf-Giant ratio (RDGR), which shows an increase by 
a factor of 2.2 ± 0.4 between z = 0.54 and z = 0.13. This 
is equivalent to an increase of 1.46 ± 0.14 in the relative 
y-band luminosity (or stellar mass) in faint red galax- 
ies with My < —20 compared to brighter systems over 
this period. This increase means that in local clusters, 
the luminosity contributed by giant and dwarf galaxies 
is comparable, whereas at z = 0.5 the giants were the 
dominant population on the color-magnitude relation. 

Our results show that there is significant evolution 
since z ~ 0.5 in the faint passive galaxy population in 
a well-defined sample of X-ray lu minous clusters . This 
agr ees with the ear l y resu lts from lDe Lucia et al.l (2004) 
and IKodama et ail (|2004h on red galaxies in a more di- 
verse range of structu res. How e ver, t his is in disagree- 
ment with the work of lAndreonl (|2006h who sees no such 
evolution. This disagreement may be simply due to An- 
dreon's use of a single cluster, as our analysis shows a 
large cluster to cluster scatter but with large errors (due 
to an uncertain field correction). We conclude that a 
large proportion of the passive galaxy population at the 
faint end of the color-magnitude sequence in local clus- 
ters either did not reside in similar, high-density envi- 
ronments 5-Gyrs ago (at z ~ 0.5) or if they were present 
in these regions then they had significantly bluer colors 
(suggesting they were actively star forming) and so do 
not fall within the color-magnitude relation. 

Clusters in the mass range studied in this work are ex- 
pected to have roughly doubled their masses since z ~ 0.5 
(Tormen 1998) and hence many of the faint red galaxies 
(or at least their progenitors) may have arrived in the 
core regions of the cluster between z ~ 0.5 and z ~ 0.1. 
However, we believe that the major driver of the evolu- 
tion we see is the transformation of blue, star-forming 
galaxies into passive, red systems which lie on the color- 
magnitude relation. If correct, this suggests that there 
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will be an increasing diversity in the star formation his- 
tories of passive galaxies at My < —20 in intermediate 
rcdshift clusters (at z ~ 0.3-0.4), and studies of age- 
sensitive indicators at these depths may uncover evidence 
for r ecent star-formati on activity within these galaxies 
(e.g. iSmail et al.ll200lh . 

We end by noting that studies such as this one can 
be extended to a wider range of environments and red- 
shifts out to z ~ 1 using the dat a from the UKIRT In- 
frared Deep Sky Survey OJKIDSS. [Lawrence et al1l200l 
in concert with the XMM Newton Large Scale Structure 
Survey (XMM-LSS). Such studies will allow the evolu- 
tion of the passive population on the color-magnitude 
relation and its build-up to be tracked as a function of 
environment and epoch and will demonstrate the impor- 



tance of including environmental effects when modelling 
the color-magnitude relation in galaxy evolution models. 
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TABLE 1 

Details of the cluster samples used in our analysis. 



Cluster 






R.A. 

■ 


Dec. 

(J2000) 






z 


Lx 

(lO^ergs" 1 ) 


tied 




DGR 


MACS 2 ~ 0.5 Sample 






























MACS J0025.4- 


-1222 


00 


25 


15 


.84 


-12 


19 


41 


0. 


.478 


12.4 


168+13 


1 


.56+0. 


35 


MACS J0257.6- 


-2209 


02 


57 


07 


.96 


-23 


26 


08 


0. 


.504 


15.4 


157+13 


1 


.47+0. 


34 


MACS J0647.7H 


-7015 


06 


47 


51 


.45 


+70 


15 


01 


0. 


584 


21.7 


183+14 


1 


.39+0. 


29 


MACS J0717.5H 


-3745 


07 


17 


31 


.8.3 


+37 


45 


05 


0. 


548 


27.4 


321+18 


1 


.25+0. 


20 


MACS J0744.8-I 


-3927 


07 


14 


51 


.98 


+39 


27 


35 


0. 


686 


25.9 


173+13 


1 


.00+0. 


22 


MACS J0911.2H 


-1746 


09 


11 


10 


.23 


+17 


46 


38 


0. 


506 


13.2 


169+13 


1 


.54+0. 


34 


MACSJ1149.5-I 


-2223 


1 1 


49 


34 


.81 


+22 


24 


13 


0. 


544 


17.3 


26o±16 


1 


.43+0. 


25 


MACS J1423.84 


-2404 


14 


2.3 


47 


.95 


+24 


04 


59 


0. 


544 


15.0 


155+13 


1 


.11+0. 


25 


MACS J2129.4- 


-0741 


21 


29 


25 


.38 


-07 


41 


26 


0. 


570 


16.4 


194+14 


1 


.27+0. 


26 


MACS J2214.9- 


-1359 


22 


14 


56 


.51 


-14 


00 


1 i 


0. 


495 


17.0 


215+15 


1 


.25+0. 


21 


LARCS z ~ 0.1 Sample 






























Abell 22 




00 


20 


38 


.64 


-25 


43 


19 


0. 


.142 


5.3 


220+21 


2 


.55+0. 


61 


Abell 550 




05 


52 


51 


.84 


-21 


03 


54 


0. 


099 


7.1 


277+22 


2 


.05+0. 


13 


Abell 1084 




10 


44 


30 


.72 


-07 


05 


02 


0. 


132 


7.4 


111+18 


4 


.00+1. 


79 


Abell 1285 




11 


30 


20 


.64 


-14 


34 


30 


0. 


106 


5.45 


391+25 


2 


.94+0. 


53 


Abell 1437 




12 


00 


25 


.14 


+03 


21 


04 


0. 


134 


7.7 


376+25 


2 


.04+0. 


35 


Abell 1650 




12 


58 


11 


.76 


-01 


45 


22 


0. 


084 


7.8 


182+20 


3 


.00+0. 


.88 


Abell 1651 




12 


59 


24 


.00 


-04 


11 


20 


0. 


085 


8.3 


232+21 


4 


.52+1. 


24 


Abell 1664 




13 


03 


44 


.16 


-24 


15 


22 


0. 


.128 


5.34 


127+18 


2 


.49+0. 


91 


Abell 2055 




15 


18 


41 


.28 


+06 


12 


40 


0. 


102 


4.8 


201+21 


4 


.54+1. 


37 


Abell 3888 




22 


34 


32 


.88 


-37 


43 


59 


0. 


153 


14.5 


124+19 


1 


.89+0. 


68 


Additional Clust- 


srs 






























CI J0152-1357 




01 


52 


43 


.91 


-13 


57 


21 


0. 


.831 


5.0 


276+17 





.77+0. 


14 


MACS J0451.9H 


-0006 


04 


51 


54 


.63 


+00 


06 


18 


0. 


430 


10.4 


179*±14 


1 


.41+0. 


30 


MACS J0712.34 


-5931 


07 


12 


20 


.45 


+59 


32 


20 


0. 


328 


6.8 


90*±10 


1 


.63+0. 


44 


CI J1226. 9+3332 


12 


26 


58 


.13 


+33 


32 


49 


0. 


890 


20.0 


232+15 


1 


.16+0. 


22 


Abell 1703 




13 


15 


00 


.70 


+51 


49 


10 


0. 


258 


8.7 


94t±10 


3 


.48+1. 


13 


MACS J1354^ 


-7715 


13 


54 


19 


.71 


+77 


15 


26 


0. 


397 


8.2 


156t±13 


1 


.66+0. 


37 



1) N ret j is the number of galaxies on the red sequence down to My = —17.75 within a 600 kpc radius of cluster centre except: * within a 
400 kpc radius, f within a 300 kpc, due to the size of the ACS image. The limited radius of M ACS 1354 is due to size of A CS image and 
flaring on the image from a bright star. The LARCS redshifts are from Pimbblct ct al. (2006). MACS redshifts come from lEbering et al.l 
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TABLE 2 

The best-fitting parameters for the luminosity function of red sequence galaxies for the MACS and LARCS clusters. 





x7" 


a 


Schechter 


4>* 


(M v ) 


Gaussian 1 

17 M 


Amp. 


LARCS 


0.28 


-1.11 ±0.05 


-21.10 ±0.11 


111. 1 ±22.2 










0.38 


-0.92 ±0.14 


-19.43 ±0.17 


331.6 ±57.2 


-20.39 ±0.12 


0.89 ±0.05 


170.5 ±21.8 


MACS 


1.29 


-0.91 ±0.02 


-21.39 ±0.05 


215.1 ± 10.9 










1.14 


-0.94 ±0.19 


-20.08 ±0.67 


223.0 ±75.4 


-20.86 ±0.16 


0.83 ± 0.08 


177.2 ±33.9 



1) Where the Gaussian fit parameters are not given this is a purely Schechter function fit to the luminosity function. 




Fig. 1. — Upper ten panels: The individual (V555 — ^814) color-magnitude diagrams for the MACS clusters. Bottom left panel: The 
combined color-magnitude diagram for the MACS sample (corresponding to restframe (U — V)—V), all clusters have been K-corrected to 
2 = 0.54. The solid lines show the limits used to define and select the red sequence in the combined clusters and the dotted line is the 5-cr 
limit / = 25.4. Bottom right panel: The combined (B — R)-R (restframe (U — V)—V) color-magnitude diagram for the LARCS sample. 
The dotted line denotes the 5-<r limit of R = 20.92 and all clusters have been K-corrected to z = 0.13. Again the solid lines show the limits 
used to select the red sequence in the combined clusters, these are transformed from the equivalent bou ndaries for the M ACS sample as 
described in the text. Color-magnitude plots for individual clusters in the LARCS survey can be found in Pimbblct et al. (2002). 
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Fig. 2. — left panel: The luminosity functions in the restframe V-band for the red sequence galaxies in the combined MACS and LARCS 
samples (normalized to the bright-end of the LARCS sample). We plot the background-corrected data and our best-fit Gaussian+Schechter 
function fits. The errors are a combination of the Poisson uncertainty and the field correction error. Note the excess of faint red galaxies in 
the lower redshift LARCS sample, compared to the more distant MACS sample, right panel: The variation in the red sequence Dwarf-Giant 
Ratio (RDGR) with redshift for clusters in our two samples within 600 kpc of the cluster centre and brighter than My = —17.75. We also 
plot the weighted mean values for each of the MACS and LARCS sample. A fit of the form (1 + z)~@ to the MACS and LARCS points 
is plotted and yields (3 = 2.5 ± 0.5. For comparison we show the equivalent measures for six additional high redshift clusters and Coma 
(RDGR= 2.8 ± 1), which follow the same trend. 



